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ABSTRACT 



Context. This paper presents the X-ray properties of a flux- and volume-limited complete sample of 16 Giga-Hertz 
Peaked Spectrum (GPS) galaxies. 

Aims. This study addresses three basic questions in our understanding of the nature and evolution of GPS sources: a) 
What is the physical origin of the X-ray emission in GPS galaxies? b) What physical system is associated with the 
i-H ' X-ray obscuration? c) What is the "endpoint" of the evolution of compact radio sources? 

Methods. We discuss the results of the X-ray spectral analysis, and compare the X-ray properties of the sample sources 
| with radio observables. 

Results. We obtain a 100% (94%) detection fraction in the 0.5-2 keV (0.5-10 keV) energy band. GPS galaxy X-ray 
spectra are typically highly obscured {(N§ PS ) — 3 x 10 22 cm -2 ; <tn h — 0.5 dex). The X-ray column density is higher 
than the HI column density measured in the radio by a factor of 10 to 100. GPS galaxies lie well on the extrapolation to 
high radio powers of the correlation between radio and X-ray luminosity known in low-luminosity FR I radio galaxies. 
On the other hand, GPS galaxies exhibit a comparable X-ray luminosity to FR II radio galaxies, notwithstanding their 
much higher radio luminosity. 

Conclusions. The X-ray to radio luminosity ratio distribution in our sample is consistent with the bulk of the high- 
energy emission being produced by the accretion disk, as well as with dynamical models of GPS evolution where X-rays 
. are produced by Compton upscattering of ambient photons. Further support for the former scenario comes from the 

location of GPS galaxies in the X-ray to 0[m] luminosity ratio versus Nh plane. We propose that GPS galaxies are 
young radio sources, which would reach their full maturity as classical FR II radio galaxies. However, column densities 
JjlO 22 cm -2 could lead to a significant underestimate of dynamical age determinations based on the hotspot recession 
velocity measurements. 



Key words. Galaxies: jets - Galaxies: active - X-ray: galaxies 

1. The nature of GPS radio galaxies tive nucleus (jElvis et al. 1994a| . Deep Chandra and XMM- 

Newton observations of GPS galaxies are scanty. One of the 
This paper presents an X-ray study of a complete radio- f exceptions is a deep XMM-Newton pointing of 3C301.1 
selected sample of giga-hertz peaked spectrum (GPS) flQ'Pea et al. 200611 : it revealed a hard X-rav emission com- 
galaxies. GPS sources are characterised by a simple con- ponent; which could be associated with hot gas shocked by 
vex radio spectrum peaking near 1 GHz (Stanghelhm 2006, the expansion of the radio source or to synchrotron self- 
Lister 2003, O'Dea 1998). They represent about 10% of Compton emission. Analysis of small samples of GPS galax- 
the 5-GHz selected sources. About half of the known i es observed with XMM-Newton were presented by Vink et 
GPS sources are morphologically classified as galaxies, a l. (2005) and Guainazzi et al. (2006). Our paper represents 
the remaining as quasars. They often exhibit symmetric, an ex t e nsion of their results 
very compact (10-100 pc) structures, reminiscent of those 

present in extended radio galaxies on much larger scales. Understanding the origin of high-energy emission in 
Little is known about their high-energy emission. GPS these objects may have important implications for the birth 
galaxies are rather elusive in X-rays (O'Dea et al. 1996). and evolution of the "radio power" in the Universe. GPS 
X-ray spectroscopic studies prior to modern X-ray observa- sources were originally suggested to represent radio galax- 
tories were inconclusive on whether this low detection rate ies in the early stage of their life (typical ages < 10 4 years; 
is due to intrinsic weakness or to obscuration of the ac- Fanti et al. 1995; Murgia 2003). This possibility was re- 
cently supported by the detection of hotspot proper mo- 
Send offprint requests to: M.Guainazzi tions (Poladitis & Conway 2003; Gugliucci et al. 2005). 
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Alternatively, as originally suggested by Gopal-Krishna & 
Wiita (1991), GPS sources could remain compact dur- 
ing their whole radiative lifetime, because interaction with 
dense circumnuclear matter impedes their full growth. 

In order to address the above issues, and provide the 
best possible estimate of the gas density in the GPS galaxy 
nuclear environment, we have undertaken an XMM-Newton 
observation program of a radio-selected complete sample 
of GPS galaxies. The three main issues which originally 
motivated our study, that will be discussed throughout the 
paper, are: 

— What is the physical origin of the X-ray emission in GPS 
galaxies? 

— Which physical system is associated with the X-ray ob- 
scuration? 

— What is the "endpoint" of the evolution of compact ra- 
dio sources? 

This paper is structured as follows: in Sect. 2 we present 
the sample. The data reduction procedure and the results 
of the spectral analysis on the unpublished sources are pre- 
sented in Sect. 3 and 4, respectively. This paper does not 
discuss the X-ray spectral ana l ysis of sources pu blished else- 
where (|Guainazzi et al. 20061 lYink et al. 2005)) . In Sect. 5 
we summarise the X-ray properties of the whole GPS sam- 
ple, and compare these with a control sample of large-scale 
radio galaxies. Finally, we discuss our results in Sect. 6, and 
summarise our findings and their possible implications for 
the evolution of GPS sources in Sect. 7. We adopt through- 
out the paper a cosmology with Hq = 70 km s _1 Mpc -1 , 
ft A =0.73, and fi m =0.27 ( |Spergel et al. 2007[ ). 

2. The sample 

The sources discussed in this paper constitute a flux- and 
volume-limited sub-sample extracted from the complete 
radio-selected sample of GPS galaxies of Stanghellini et al. 
(1998). We refer to this sub-sample as "our GPS sample" 
hereafter. We selected all sources with redshift z<l, and 
flux density at 5 GHz >1 Jy. The whole sample (16 sources) 
has been observed with XMM-Newton across different ob- 
serving cycles. The only exceptions are PKS 0941—08 and 
PKS1345+125, for which archival Chandra and ASCA data 
are available, respectively. Preliminary results, based on a 
small sub-sample of 5 objects, were presented in Guainazzi 
et al. (20060 Some of the sample sources were included in 
Vink et al. (2005). The whole sample is listed in Table [Q A 
more general discussion of the X-ray and multiwavclcngth 
properties of our whole GPS sample is deferred to Sect. 5. 
A summary of the X-ray and radio properties of the whole 
sample is in Appendix B. 

3. Observations and data reduction 

In this paper [as also originally done by Guainazzi et al. 
(2006) and Vink et al. (2005)] we will consider only X- 
ray data taken with the XMM-Newton EPIC cameras (pn, 
IStriider et al. 20011 MOS. ITTIrner et al. 2001) . because the 

1 Guainazzi et al. (2006) present also data of 

COINSJ0029+3456; this source was later discovered to 

host a blazar, and will not be considered in the sample 
discussed in this paper. 



Table 2. Properties of the X-ray observations discussed in this 
paper. Exposure times (T exp ) and Count Rates (CR) refer to the 
pn in the 0.5-10 keV band (unless otherwise specified), r is the 
count rate threshold applied in the determination of the Good 
Time Interval for scientific product extraction (details in text), p 
is the radius of the scientific products extraction region. Nh,Go,i 
is the column density due to gas in the Milky Way along the 
line-of-sight to the GPS galaxy (Kalberla et al. 2005) in units 
of 10 20 cm" 2 . 



Source 








P 


CR 






(ks) 


(s- 1 ) 


(") 


(ht 3 s- 1 ) 


4C+00.02 


2.7 


22.2 


2.1 




a 


PKS0428+20 


19.6 


15.4 


0.75 


23 


7.1 ±0.8 


4C+14.41 


2.0 


16.3 


0.5 


15 


4.8 ±0.8 


4C+32.44 


1.2 


20.5 


0.5 


16 


25.6 ± 1.2 


PKS1607+26 


3.8 


18.5 


18 


20 


39 ±7 


PKS2008-068 


5.0 


23.2 


0.75 


15 


2.0 ±0.4 


PKS2127+04 


5.0 


10.8 


0.5 


25 


21.2 ± 1.9 



"detected only in the soft X-ray band (0.5-2 keV), with a count 
rate: (2.3 ± 0.7) x 10" 3 s _1 . 



sources were too faint to yield a measurable signal in the 
high-resolution RGS cameras. Observational information 
for the sources for which new measurements are presented 
here is listed in Table O 

XMM-Newton data were reduced with SASv7.1 
(|Gahriel et al. 2003)) according to standard procedures as 
in, e.g., Guainazzi et al. (2006). The most updated cali- 
bration files available at the date of the analysis (February 
2008) were used. Source scientific products were accumu- 
lated from circular regions surrounding the position of the 
optical nucleus of each source (extracted from the NED 
catalogue j3- The size of the source extraction regions are 
shown in Tabled] Following Guainazzi (2008), background 
scientific products were extracted from source free circu- 
lar regions close to the source and on the same CCD as 
the source for the MOS cameras; and from source-free re- 
gions centred at the same row in detector coordinates as the 
source in nearby CCDs for the pn. As many of the sources 
were X-ray faint, particular care has been applied in the 
choice of flaring particle background rejection thresholds 
optimising the signal-to-noise ratio of the final scientific 
products. Using a single-event, E >10 keV full field-of-view 
light curve as a monitoring tool of the instantaneous inten- 
sity of the background, ten different logarithmically spaced 
thresholds between 0.1 s _1 and about two times the high- 
est light curve count rates were tried. For each threshold 
the radius of the source extraction region was also varied to 
obtain the highest number of net counts for a given signal- 
to-noise ratio. Spectra were binned in such a way as to 
avoid oversampling of the intrinsic instrumental energy res- 
olution by a factor larger than 3, and to have at least 25 
background-subtracted counts in each spectral bin. These 
conditions ensure the applicability of the x 2 statistic as a 
goodness-of-fit test. 

In this paper, errors on the spectral parameters and on 
any derived quantities are at the 90% confidence level for 
one interesting parameter; errors on the count rates and 
derived quantities are at la level. Whenever statistical mo- 
ments or correlations on distributions including upper lim- 

2 http : //nedwww. ipac . caltech.edu/ 
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Table 1. The GPS sample discussed in this paper. 
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NED name 


J2000 name 




z a 


Observation ID b 


Observation date 


Reference 2 


4C+00.02 


002225+001456 





.305 


0407030101 


2006-Jun-26 


(1) 


COINSJ0111+3906 


011137+390628 





.668 


0202520101 


2004-Jan-09 


(2) 


PKS0428+20 


043103+203734 





.219 


0407030201 


2006-Sep-02 


(1) 


PKS0500+019 


050321+020305 





.585 


0205180601 


2004- Aug- 18 


(3) 


B30710+439 


071338+434917 





.518 


0202520201 


2004-Mar-23 


(2) 


PKS0941-080 


094336-081931 





.228 


Chandra 




(3) 


B1031+567 


103507+562847 





.450 


0202520301 


2004-Oct-21 


(2) 


4C+14.41 


112027+142054 





.362 


0502510201 


2007-Jun-13 


(1) 


4C+32.44 


132616+315409 





.370 


0502510301 


2007-Dec-05 


(1) 


PKS1345+125 


134733+121724 





.122 


ASCA 




(4) 


4C+62.22 


140028+621038 





.431 


0202520401 


2004- Apr- 14 


(2) 


OQ+208 


140700+133438 





.077 


0140960101 


2003-Jan-31 


(•>) 


PKS1607+26 


160913+264129 





.473 


0502510401 


2008-Jan-17 


(1) 










0502510801 


2008- Jan- 19 


PKS2008-068 


201114-064403 





.547 


0502510501 


2007-Oct-12 


(1) 


PKS2127+04 


213032+050217 





.990 


0502510601 


2007-May-17 


(1,6,7) 


COINSJ2355+4950 


235509+494008 





.238 


0202520501 


2003-Dec-26 


(2) 



a Data taken from NED 

b XMM-Newton observatio n ID 

c (1) this paper; (2) IVink et al. 20051 (3) IGuainazzi et al. 20061 (4) IQ'Dea et al. 20001 (5) IGuainazzi et al. 20041 
(6) ISiemiginowska et al."2005| (7)[S iemiginowska et al. 2008 



its are calculated, an extension of the regression method 
on censored data originally described by Schmitt (1985) 
and Isobe et al. (1986) has been used. More details on this 
method were presented by Guainazzi et al. (2006). 

4. Results 

In this Section we present spectral-analysis for the 7 un- 
published sources in our GPS sample. 

No significant variability in cither integrated X-ray flux 
or spectral shape was detected in any source presented in 
this paper on timescales ^10 4 s. We therefore focus on the 
properties of their time-averaged spectra. 

For 3 of these sources the number of degrees of free- 
dom in the binned spectra was larger than 4: 4C+32.44, 
PKS1607+26, PKS 2127+04. In these cases a standard 
spectral analysis was possible. The spectra were fitted in 
the 0.2-10 keV energy range with Xspec (version 11; 
lArnaud r9 96 ). A model consisting of three components was 
used: 

e -N H}Gal a(E) x e -N„a[E(l+z)] x K £-T Q\ 

where the photo-electric absorption components use 
Wisconsin cross-sections (Morrison & McCammon 1983), 
and K is the unabsorbed spectral normalisation at 1 keV. 
We'll refer to this model as our "baseline" model hereafter. 
The column density Nn,Gai was kept fixed to the contri- 
bution by intervening gas in the Galaxy as measured in 
the Leiden/ Argentine/Bonn (LAB) Survey of Galactic HI 
(Kalberla et al. 2005[rl , whereas Nh is the column density 
of local gas at the galaxy redshift. The results of the spectral 
fitting are shown in Fig. [1] and summarised in Table [31 X- 
ray luminosities in this table and hereafter are K-corrected 
and based on the best-fit model for each source. Addition 
of a Gaussian profile to the baseline model to characterise 

3 http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 



line emission never improved significantly the quality of the 
fit. The upper limits on the Equivalent Width (EW) of an 
unresolved Fe K a neutral fluorescent line are also reported 
in Table[3j They are generally inconclusive. Best-fit models 
and residuals are shown in Fig. [TJ 

For 4 other sources (4C+00.02, PKS 0428+20, 
4C+14.41, and PKS 2008-068) the low signal-to-noise did 
not allow for spectral analysis. We have therefore based our 
estimation of the spectral parameters of the baseline model 
on the Hardness Ratio (HR), here defined as the ratio be- 
tween the counts in the energy bands 1-10 keV and 0.2- 
1 keV. The measured HRs (or lower limit thereof) have been 
compared with the predictions of grids of simulated baseline 
models. Iso-HR contour plots in the T ([0.5:3]) versus Nh 
([10 20 ,10 24 cm -2 ]) parameter space were built (see Fig. [5]) 
The confidence interval in column density was estimated 
as the minimum and maximum value of the 1.6cr interval 
around the curve corresponding to the nominal HR, when T 
was constrained in the range: [0.63,2.62]. The photon-index 
range corresponds to the (r) ± 3<r interval of the T dis- 
tribution, calculated on the whole sample of GPS galaxy 
for which the data quality allowed us to perform a spec- 
tral analysis (cf. Fig. [6]). The resulting constraints on the 
column density are shown in Table [H 

4.1. A Compton-thick AGN in PKS1607+26? 

The fit of the EPIC spectra of PKS1607+26 yields an un- 
usually flat spectral index: T = 0.4 ± 0.3. A flat spec- 
trum may be indicative of a blazar-type spectral compo- 
nent. Alternatively, in radio-quiet AGN spectral indices 
r ~ are generally interpreted as evidence for a Compton- 
thick AGN, whose primary X-ray emission is totally sup- 
pressed by optically thick matter with a column density 
Nh > cr^ 1 ^ 1.6 x 10 24 cm" 2 (see Comastri 2004 for a re- 
view). In Compton-thick AGN, residual X-ray emission red- 
wards the photoelectric cutoff could be due to Compton- 
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Table 3. Best-fit parameters and results for the sources where spectral analysis was possible. 



Source 


N H 
(10 21 cm 2 ) 


r 


F2-10 

(10 -13 ergs cm" 2 s" 1 ) 


log(LI-io) 


Fe K a EW 
(keV) 


X'/v 


4C+32.44 


i.2±8:S 


1.74 ±0.2 


0.78 ±0.13 


43.57±S:S 


<1.0 


20.6/36 


PKS1607+26 


<2 


0.4 ±0.3 


4.2±J:§ 


44.27 ±0.01 


<0.7 


32.2/21 


PKS2127+04 


< 19 


1.98 t° ; 5 4 


0.49±[J;i? 


44.46 ± 0.06 


<0.9 


14.7/16 



a in erg s 




Fig. 1. Spectra of the three EPIC instruments (upper panels), and residuals in units of standard deviation (lower panels) 
when the baseline model is applied. In addition to the binning applied to the spectra to ensure the applicability of the 
X 2 goodness-of-fit test, the spectra were further rebinned in this Figure such that each spectral channel corresponds to 
a signal-to-noise ratio >3 for plotting purposes only. Green: pn; black: MOS 1; red: MOS 2. 



Table 4. Constraints on the intrinsic column density derived 
from iso-HR contours in the F versus Nh parameter space. No 
constraint can be derived for 4C±00.02. 



Source 


Nn,max 


-F2-10 


log(L 2 -io) 




(10 21 cm" 2 ) 




PKS0428±20 


<6.9 


[0.3, 0.6J 


43.1 


4C+14.41 


<1.6 


[0.2, 0.5] 


43.1 


PKS2008-068 


<4.8 


[0.2, 0.5] 


43.8 



"ranges in units of 10 -13 erg s _1 cm -2 corresponding to the 
Nh range as in this table, and to the photon index in the range: 
[0.63,2.62] (more details in the text). 



reflection of the otherwise invisible primary radiation off the 
obscuring matteiQ. Fe K a fluorescent emission with large 
EW 600 eV; Risaliti 2002, Guainazzi et al. 2005) is 
the unmistakable sign of a Compton-thick AGN. We have 
therefore applied the baseline spectral model for Compton- 
thick AGN to the PKS1627±06 pn spectrum (the MOS 
spectra have very poor statistics, and do not provide any 
further constraints): a pure Compton-reflection continuum 
(model pexrav in Xspec; Magdziarz & Zdziarski 1995) plus 
a Gaussian unresolved emission line. In order to take ad- 
vantage of the full energy resolution of the EPIC cameras 
we fit the unbinned spectrum using the Cash statistic, C 
(jCash 1976[) . Although the Cash statistic does not provide 
an absolute level of statistical confidence on the quality of 
the fit, the resulting value (C=25.2/21 degrees of freedom) 
and the smoothness of the residuals (Fig. [3]) indicate that 

4 the GPS galaxy OQ±208 was the first radio-loud Compton- 
thick AGN discovered (|Guainazzi et al. 20"04|) 



PKS1607+26-EPIC/pn 




Energy (keV) 

Fig. 3. EPIC-pn spectrum of PKS1607±26 in the 2-10 keV 
energy band (observer's frame). The data points have been 
rebinned such that each displayed spectral channel has a 
signal-to-noise ratio larger than 1 . The solid line represents 
the best-fit reflection-dominated model (details in text). 



the fit is good. Once this continuum is adopted, an emission 
line is required at a confidence level larger than 90% for two 
interesting parameters (or 95% for one parameter; Fig. 0]). 
The EW (500 ± 300 eV) is consistent with values typically 
observed in Compton-thick AGN. However, the best-fit cen- 
troid energy of this feature (E* =0 = 4.12±g;g| keV) is 
slightly inconsistent with Fe K Q fluorescence if converted 
into the rest frame with the quoted redshift for this object: 
E* =0A73 — 6.06±Qgy keV. No other known atomic transi- 
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4C+00.02 - XMM/EPIC - HR> 0.91 PKS0428+20 - XMM/EPIC - HR— 9.05±3.96 




log(N H , cm" 2 ) log(N H , cm" 2 ) 

Fig. 2. Confidence level loci in the (r,7V#)-plane. For 4C+00.02 the white area shows the locus of the parameter space 
where the model is consistent with the measured HR. In the other panels the grey scale indicates the deviation from the 
measured HR in units of standard deviations. 



, PKS1607+26-EPIC/pn 

6 i 1 1 




Q.I 
c 

.S i- - 

- 1 o I i , i , i ; , 

x 3.8 4 4.2 4.4 

Line Energy (keV) 

Fig. 4. Iso-x 2 contours in the centroid energy (observer's 
frame) versus normalisation plane for the Gaussian profile 
in the best-fit reflection-dominated model for PKS1607+26. 
The contours correspond to the 68%, 90% and 99% confi- 
dence level for two interesting parameters, respectively. 



Philips & Shaffer (1983). Biretta et al. (1985) refer to pre- 
vious mis-identification of this object, and indeed the value 
of z=0.47 in the original paper corresponds to the wrong 
identifier in Biretta et al. (1985). If the emission line is in- 
deed the Fe K Q , the redshift is z ~0.55. Notwithstanding 
the ultimate redshift value, the very flat hard X-ray spec- 
tral index, alongside the detection of a strong emission-line 
features point to a Compton-thick nature for this object. 
In order to determine a lower-limit to the column density 
covering the primary continuum, we have added to the pure 
reflection model a power law obscured by photoelectric ab- 
sorption. We have used a combination of the models phabs 
and cabs in Xspec to take into account the contribution 
of Compton scattering, which is no longer negligible for 
column densities close to the Compton-thick limit. The ob- 
scuring column density is constrained to be higher than 
6 x 10 23 cm -2 (90% confidence level for one interesting 
parameter; Fig. [5]). This is the value we will associate with 
PKS1607+26 in all subsequents plots and distribution func- 
tions. However, given the pending uncertainties on the ul- 
timate identification of this source, we will not use this 
measurement to derive statistical estimators in the Nh dis- 
tribution of the whole sample. 



tion could be responsible for a feature at this energy We do 
not have a convincing explanation for this finding. The red- 
shift quoted in the literature for PKS1607+26 dates back to 
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o PKS1607+26-EPIC/pn 




Fig. 5. Iso-x 2 contours for the column density covering a 
power-law primary continuum derived from the EPIC pn 
spectrum of PKS1607+26 in the reflection-dominated sce- 
nario. The lines correspond to the 68%, 90% confidence 
level for one interesting parameter and to the 90% confi- 
dence level for two interesting parameters, respectively. 



5. Comparison with control samples of "normal" 
radio galaxies 

In this Section we present X-ray spectral properties of the 
flux density and redshift-selected complete sub-sample ex- 
tracted from the Stanghellini et al. (1998) GPS sample de- 
scribed in Sect. 2. Readers are referred to the Guainazzi et 
al. (2006) and Vink et al. (2005) for the spectral analysis of 
the sources not discussed in this paper. We have repeated 
the analysis of the sources of the Vink et al. (2005) with the 
same reduction and data screening criteria as in Guainazzi 
et al. (2006) and in this paper. The results of our re-analysis 
are consistent with theirs. 5 GHz luminosities are taken 
from Stanghellini et al. (1998) and O'Dea (1998). 

Our goal is also to compare the properties of the 
complete radio-selected flux-limited GPS sample with a 
control sample of "normal" radio galaxies. The control 
sample has been built from results recently published 
in the literature. It is based on observations of z < 1 
radio galaxies taken by ASCA (jSambruna et al. 1999[) . 
BeppoSAX IjGrandi et al. 2006P, XMM-Newton and 
Chandra (FR I: IDohato et al. 20041 lEvans et al. 20061 
IBalmaverde et al. 20061 FR II IBelsole et al. 20061 
lEvans et al. 20061 IHardcastle et al. 2006|) . Only one 
measurement per source has been retained in the control 
sample, based on the latest published result. However, we 
have considered the latest Chandra measurement, even 
when a later XMM-Newton observation was available, 
under the assumption that the superior spatial resolution 
of the Chandra optics provides a more reliable measure- 
ment of the core emission. The control sample comprises 
93 sources (32 FR I, 54 FR II, the remaining ones have 
no FR classification). We stress that the control sample 
is neither complete nor unbiased. Moreover, it is not well 
matched with our GPS sample in redshift. The probability 
that the redshift distribution of the GPS sample is the 
same as in the entire control sample is 2% only. This 
difference is mainly due to FR Is being generally at lower 
redshift than our GPS sample, whereas our GPS and the 
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1 - 

^ 
0.5 




1.0 



1.5 2.0 
Photon Index 



2.5 



3.0 



Fig. 6. Distribution of the photon index for the GPS sub- 
sample (8 objects), where data quality warranted spectral 
analysis. 



control FR II samples are well matched in redshift (cf. 
Fig. [5]). Whenever pertinent, we will explicitly show in 
the following the redshift dependence of the observables, 
and discuss any possible bias associated with comparing 
samples inhomogeneous in redshift coverage. 



5.1. X-ray detection fraction 

We obtain a very large detection fraction; all the sources 
of our sample yield a detection in the soft X-ray band 
(0.5-2 keV), whereas 15 out of 16 are detected in the full 
band (0.5-10 keV). All of them but one (PKS1345+125) 
were unknown in X-rays prior to our Chandra (see also 
Siemiginowska et al. 2008) and XMM-Newton observations. 



5.2. Spectral shape 

In Fig. [6] we show the distribution of spectral indices for 
the 8 GPS of our sampl^l, where the number of counts is 
good enough for the spectral analysis to be possible. The 
distribution has a mean value (r) = 1.66 ± 0.40, and a 
standard deviation or = 0.36. The weighted mean is T = 
1.61 ± 0.05 if the measurements are weighted according to 
the inverse square of their statistical uncertainties. 



5.3. Obscuration 

In Fig. [7] the distribution of column density for the GPS 
sources of our sample is shown. The mean value of the dis- 
tribution is (N§ PS ) = 3 x 10 22 cm with a standard devia- 
tion cr/Vff— 0.5 dex. The same Figure shows the comparison 
with the control sample. The core emission of FR I radio 



5 4C+32.44, 4C+62.22, B03710+439, COINSJ2355+4950, 
PKS0500+019, PKS1345+125, PKS2127+04, OQ+208 



Tengstrand et al.: X-ray view of GPS radio galaxies 



7 




12 3 4 

log (N H , 10 2 °cm 2 ) 



Fig. 7. Distribution of the core obscuring column densities 
for the GPS sample (bottom panel), the FR I (upper panel) 
and FR II (middle panel) control samples. Shaded areas 
indicate upper limits; empty areas indicate lower limits. 
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Fig. 8. Core obscuring column density as a function of 
redshift for the GPS galaxies (filled dots) and the FR II 
control sample (empty squares). 



galaxies is generally unobscured or only mildly obscured. In 
the Donato et al. (2004) FR I sample, less than one-third of 
the sample exhibits excess obscuration above the Galactic 
contribution, with rest-frame column densities in the range 
10 20-21 cm -2 , thus significantly lower than observed in our 
sample. The average of the column density distribution in 
the FR I control sample is (iV^ RI ) = (3.3±^) x 10 21 cm~ 2 . 

FR II cores tend instead to include a heavily obscured 
component. However, a detailed quantitative comparison 
between the GPS and the FR II control sample is difficult. 
For 11 out of 54 FR II sources neither a measurement nor 
an upper limit on the column density is available in the 
literature. These sources are generally considered as unob- 
scured (Hardcastlc et al. 2006). Lack of inclusion of these 
sources could potentially bias the control sample Njj distri- 
bution towards higher values. Bearing this caveat in mind, 
GPS galaxies seems to fill a gap in the Nh distribution 
between highly obscured (iJlO 23 cm -2 ) and unobscured 
(^10 22 cm~ 2 ) FR II spectral components. 

A potential area of concern is the comparison of col- 
umn density measurements in samples, which are not well 
matched in redshift. However, Fig. [H] shows that this bias 
is not responsible for the difference between the average 
of the column density distribution in the GPS and FR II 
samples. Moreover, low-redshift GPS galaxies exhibit col- 
umn densities not systematically lower than high-redshift 
ones. 

Hardcastle et al. (2006) remark that heavily absorbed 
nuclei are rather common in narrow-line radio galaxies, 
whereas they are comparatively rare in Low-Excitation 
Radio Galaxies (LERG; |Laing et al. 1994fl . There are 7 
LERGs in our control sample; 5 of them have no column 
density measurement; the remaining two have column den- 
sity of ~3 x 10 22 cm' 2 (3C 123; IHardcastle et al. 2006P 



and ~10 23 cm" 2 (3C 427.1; IBelsole et al. 2006ft . Taking 
into account the low number statistics and the uncertain- 
ties on the column density upper limits on formally "unob- 
scured" LERGs, the comparison between X-ray obscuration 
in LERGs and GPSs is inconclusive. 

The X-ray column density is significantly larger than the 
column density measured by HI observations by a factor 10 
to 100. The comparison is shown in Fig. [Sj The estimate 
of the HI column density depends on the values assumed 
for the spin temperature, T s , and for the fraction of back- 
ground source covered by the absorber, f c . The data in 
Fig.[Uassume T s =100 K and / c =l (jVermeulen et al. 2 003. 



IPihlstrom et al. 20031 |Gupta et al. 2006| ). The X-ray ver- 
sus HI column density relation can be fit with a zero-offset 
linear function if T s is of the order of a few thousands K 
(jOstorero et al. 2009|) . Alternatively, a low covering frac- 
tion could be responsible for the large X-ray to HI column 
density ratio, although this explanation is less likely given 
the large column densities measured in X-rays. 

Holt et al. (2003) proposed an "onion-skin" model for 
the nuclear environment gas in B1345+125, to explain the 
reddening properties of the different components of the op- 
tical lines. The jet would pierce its way through a dense 
cocoon of gas and dust, with decreasing density at larger 
distances from the radio core. If this scenario would apply 
to the whole class of GPS samples, and X-ray and radio 
emission originate in the same physical system, one might 
expect an anti-correlation between the measured column 
density and the size of the radio structure, with a large 
scatter due to the unknown line-of-sight angles distribu- 
tion. This correlation is shown in Fig. [TDJ A censored fit 
on this data with a function: \og(Nn) = A + B x \og(rk pc ) 
yields: A = 21.4 ± 0.4, and B = -1.2 ± 0.3, where r kpc 
is the size of the radio structure in units of kpc. This re- 
sult is formally robust, but admittedly mainly driven by 
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Fig. 9. Comparison between the column densities mea- 
sured in X- rays (this paper) and w ith atomic hydrogen ob- 
servations (Pihlstrom et al. 2003)) . The lines represent loci 
of constant X-ray versus HI column density ratio for ratio 
values of 1, 10, and 100 (from right to left), respectively. 
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Fig. 10. X-ray column density versus the size of the radio 
structure. The solid line represent the best censored data 
linear fit with a function: \og(Nn) = A + B x log(rfc pc ); the 
dotted lines correspond to ±lrr uncertainties on the best-fit 
parameters 



the two extreme data points. A confirmation of this cor- 
relation by increasing the number of small (r <100 pc) 
and large (r > 1 kpc) objects for which spectroscopic X- 
ray data are available is a task we are actively pursuing. 
Interestingly enough, an anti-correlation between the lin- 
ear dimension of the sub-galactic radio galaxy and the ra- 
dio HI column density was discovered by Pihilstom et al. 
(2003). As already pointed out by Gupta et al. (2006), this 
anti-correlation could be driven by small sources probing 
gas closer to the AGN and hence at a higher spin tem- 
peratures. In this context, it is also interesting to observe 
that GPS quas ars exhibit no absorption (w ith upper limits 
r~ 10 21 cm~ 2 ; Sicmig mowska et al. 2008[ ), as well as dif- 



fuse emission associated with jets, binary structures or em- 
bedding clusters. The detection rate as well as the column 
density of HI absorption increases with core prominence 
( |Gupta k, Saikia"2 006). The core prominence is a statisti- 
cal indicator of the orientation of the jet axis to the line 
of sight. On average HI absorbers are more common and 
exhibit higher column densities in galaxies than in quasars. 
This can be explained if the HI absorbing gas is distributed 
in a circumnuclcar disk much smaller than the size of the 
radio emitting region, and only a small fraction of it is ob- 
scured in objects at large inclinations. 

5.4. Radio-to-X-ray correlations 

In the left panel of Fig. [TT] we compare the logarithmic ra- 
tio between the 2-10 keV and the core 5 GHz luminosity 
(L 5G Hz = v$ghz x S 5GHz , where S 5GHz is the luminos- 
ity density at 5 GHz) for the GPS and the control sample. 
Values for the GPS sample range between -0.5 and 1.5. 
No clear dependence on redshift is observed. GPS galaxies 



are X-ray under- luminous by about an order of magnitude 
with respect to their radio power once compared to FR II 
radio galaxies, blazars (Fossati et al. 1998)) and radio-loud 
quasars (Elv is et al. 1994b)) . On the other hand, the X-ray- 
to-radio luminosity ratios in GPS galaxies well match values 
observed in FR I galaxies: log(Lx / LhGHz) < 1.8. It is im- 
portant to stress again that there is almost no overlap in 
redshift between the GPS and the FR I samples, though. 

The interpretation of the X-ray-to-radio luminosity cor- 
relation depends on the origin of the bulk of the VLA radio 
emission in compact galaxies. VLBI observations of GPS 
galaxies revealed a fraction of Compact Symmetric Objects 
(CSO ) between -30% and 100% (Stang hellini et al. 1997, 
1999; |Xiang et al. 2005] ILiu et al. 2007)) . Three objects in 
our sample show a CSO morphology: PKS 0050+019, 
PKS 1345+125, and PKS 2008-068 (Stanghellini et al. 
1997, 1999), although in all these cases the morphology 
is rather complex, with multiple components on scales 
^20 pc. High-resolution, multi-frequency observations of 
our sample would be required to estimate which fraction of 
the VLA flux can be safely attributed to a core. 

From Fig.[TT]a possible bimodality of the X-ray-to-radio 
luminosity ratio in the GPS sample is apparent. The fit of 
the cumulative distribution function of this quantity with 
a single Gaussian yields a Kolmogornov/Smirnov value of 
0.33, corresponding to null hypothesis probability of about 
4%. A fit with a double Gaussian yields instead a value of 
0.18, with a null hypothesis probability of 65%. We consider 
this as a tentative piece of evidence for bimodality only. We 
will refer in the following to "X-ray bright" and to "X-ray 
weak" GPS galaxies as those, whose log(Lx)/{L 5 ghz) ra- 
tio is larger/smaller than 0.5, respectively. No significant 
difference in the spectral shape between the the "X-ray 



Tengstrand et al.: X-ray view of GPS radio galaxies 



9 




Fig. 11. Left panel: 2-10 keV versus 5 GHz logarithmic ratio versus redshift for the GPS (filled dots) and the FR II 
sample (empty squares). The obliquely shaded box indicates the locus of the FR I Chandra sample; the horizontally shaded 
box the locus of the blazar sample of Fossati et al. (1998); the dot-dashed line the locus corresponding to a typical Spectral 
Energy Distribution of a radio-loud quasar after Elvis et al. (1994b). Right panel: 2-10 keV versus 5 GHz core luminosity 
for the GPS galaxies (filled circles), and a control sample of radio galaxies (FR \:empty circles, FR II: empty squares). 
The lines represent the best-fit regression line for censored data for X-ray weak GPS (continuous), FR I (dashed) and 
FR II (dashed- dotted) 



bright" and "X-ray weak" sources was observed. In par- 
ticular, the obscuring column density distributions are in- 
distinguishable. 

In the right panel of Fig. [Tl] the 2-10 keV luminos- 
ity is plotted against the 5 GHz luminosity. In FR I 
galaxies a strong correlation between core X-ray, ra- 
dio and optical flux is known (Chiabe rge et al. 1999 
Har dcastle fc Worrall 2 000). In our control sample, the 
slope of this correlation is consistent with unity: 0.97±0.06. 
Interestingly enough, this slope is consistent with the slope 
observed in "X-ray weak" GPS galaxies 0.95 ± 0.12, with 
a ~0.5 dcx offset at face-value. The "X-ray bright" sample 
exhibits a significantly flatter slope, although with large un- 
certainties (0.1 ±0.2). Again, it is hard to estimate any bias 
associated with the coarse spatial resolution of the GPS ra- 
dio measurements. 



6. Discussion 

In this Section we will review the main X-ray observational 
properties of our whole GPS sample, trying to address the 
three issues which originally motivated our study: 



— What is the physical origin of the X-ray emission in GPS 
galaxies? 

— Which physical system is associated with the X-ray ob- 
scuration? 

— What is the "endpoint" of the evolution of compact ra- 
dio sources? 



6.1. X-ray spectral support for an accretion-disk origin 

The distribution of spectral indices in the GPS sample 
by itself does not provide any stringent constraints on 
the origin of X-ray emission in GPS galaxies. The mean 
value, (r) = 1.66 (or — 0.36), is consistent with spec- 
tral components associated with the jet in radio galaxies 
(r = 1.88±0.02), as well as with accretion (r = 1.76±0.02; 
Evans et al. 2006), although is nominally closer to the lat- 
ter. A possible clue may come from the fact that the column 
density measured in X-rays is invariably larger by 1-2 or- 
ders of magnitudes than that measured in radio. This find- 
ing could be naturally explained by X-rays being produced 
in a smaller region than the radio. On average, the radio 
morphology of compact radio sources strikingly resembles 
that of large-scale radio doubles, although on a scale which 
is entirely confined within the optical narrow-line emission 
regions. Radio emission traces therefore the radio hotspot 
and lobes. X-rays could be instead generated by the base 
of the jet, or by the accretion disk. 

Support for the X-rays arising from a relatively compact 
region comes from the comparison with radio-quiet AGN. 
Once a similar baseline model is employed, Seyfert galaxies 
have: (T Sy ) = 1.64 ±0.05 (Bianchi et al., submitted]! GPS 
galaxy X-ray spectra lack the complexity that Seyferts typ- 
ically exhibit. There is no strong evidence for a soft excess 
(with the only exception of OQ±208; Guainazzi et al. 2004), 
warm absorber, warm scattering or Fe K Q fluorescent emis- 



This is not the intrinsic spectral index, which can be signifi- 
cantly larger due to the hardening effect of ionised absorbers of 
disk/torus reflection. 
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sion (with, again, the notable exceptions of OQ+208 and, 
possibly PKS1607+26) in our sample. However, most of 
the GPS galaxy spectra collected so far with either XMM- 
Newton or Chandra do not possess the statistical quality 
that would be needed for these additional spectral features 
to be unambiguously detected. 

The scatter in the X-ray to radio luminosity ra- 
tio, and its lack of dependence on source size and 
age (cf. also Fig. [T5|) may indicate a link between ac- 
cretion disk and jet activity primarily driven by disk 
instabilities. 10-20% of GPS objects exhibit very ex- 
tended radio (|Baum et al. 19901 |Stanghellini et al. 19901 
ISchoe nmakcr s~et al. 19991 Marccki et al. 20031 or X-ray 
( jSiem iginows kaet al. 2002} 2003) emission. These com- 
ponents have been interpreted as remnants of past en- 
hanced activity. A similar behaviour on much shorter time 
scales is observed in Galactic Black Hole Candidates and 
micro-quasars, suc h as GRS 1915+105 (|Belloni et al. 2000, 
IFender et al. 2004)) . Jet blobs are supposed to be fed by 
the evacuation of the innermost accretion disk regions. This 
mechanism yields alternating X-rays- (disk-dominated) and 
radio-bright (jet blob-dominated) phases. Models based 
on disk radiation pressure instabilities reproduce well the 
timescales of these transitions (Czerny et al., in prepara- 
tion), although they do not make yet specific predictions 
on the evolution of the spectral energy distribution. 

6.2. Dilution by X-rays from radio-emitting regions? 

Most likely, the baseline model is indeed too simple. X-ray 
absorption could be "diluted" by X-rays coming from the 
high surface-brightness radio components, thus complicat- 
ing the interpretation of the results derived by our simple 
baseline model. This scenario would also explain the (still 
tentative) anti-correlation between the X-ray column den- 
sity and the size of the radio source (Fig. [T0j) . In larger 
sources, a larger fraction of the X-ray emission associated 
with the radio hot-spots or lobes may be visible beyond 
the rim of the obscuring matter. Should this indeed be the 
case, we should, however, observe deviations from a sim- 
ple power-law spectral shape. We indeed observe a soft ex- 
cess in the radio-loud Compton-thick GPS galaxy OQ+208. 
This is the closest object in our sample, and the only one 
where high-resolution spatially-resolved spectroscopy with 
Chandra could provide direct observational clues on the 
physical location of the X-ray emitting plasma. 



6.3. The X-ray obscuring environment 

Lacking other direct evidence from X-rays alone, one may 
use multiwavelength diagnostics to obtain further clues on 
the origin of the X-ray emission in GPS galaxies. A well 
established diagnostic tool for X-ray obscuration in radio- 
quiet AGN involves the comparison between the column 
density measured in X-rays and the absorption-corrected 
ratio of X-ray to 0[iii] fluxes (Mai olino et al. 1998)) . In the 
context of Seyfert galaxies, this diagram is used to calibrate 
the latter quantity as an estimator for obscuration. We 
use here the same plot with a different purpose, namely to 
test whether the ionising continuum powering the Narrow 
Line Regions in GPS galaxies has the same properties as 
in Seyfert galaxies once normalised to the X-ray primary 
emission. If this is the case, one may conclude that the 
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Fig. 12. X-ray column density versus the ratio between the 
absorption-corrected 2-10 keV and the Ofiii] fluxes. The 
filled squares are the 5 GPS galaxies in our sample for which 
O [in] measurements are available; the empty circles repre- 
sent a sample of Seyfert 2 galaxies after Guainazzi et al. 
(2005). 



"normalising primary continua" share the same properties 
in the two populations. The results of this comparison are 
shown in Fig. [T2l 5 objects in our GPS sample have spec- 
troscopic 0[iii] measurements (O'Dea 1998, Labiano et al. 
2005, 2008 and references therein). The "control sample" is 
a collection of Seyfert 2-galaxy measurements in Guainazzi 
et al. (2005). The agreement is good. Although more GPS 
data would be required to ensure a homogeneous coverage 
of this plane, this evidence further points towards an accre- 
tion origin for the X-ray emission in GPS galaxies. 

Interestingly, the fraction of Compton-thick AGN in 
the GPS sample (0.12 ± 0.09, if PKS1607+26 is con- 
sidered) is comparable to fractions observed in radio- 
quiet AGN (Hcckma n et al. 2005)) . Conversely, only one 
Compton-thi ck AGN has been detected in large-scale ra- 
dio galaxies (jErlund et al. 2008)) . although FR II typically 
exhibit spectral components with very large obscuration 
(jBelsole et al. 2006)lEvans et al. 20"06)) . 

6.4. X-ray emission entirely associated with radio-emitting 
regions? 

We now investigate the implications of the difference 
between X-ray and radio HI column densities, possibly 
due to the ionisation state of a single gas system si- 
multaneously covering the radio and the X-ray emission. 
Attributing this difference solely to ionisation effects would 
mean ionisation fractions of 90% to 99% (see the dis- 
cussion in Vink et al. 2005). Photoionisation simulations 
with CLOUDY (jFerland et al. 1998)) show that this corre- 
sponds to an ionisation parameter ^> 20 for a gaseous 
nebula photoionised by a typical AGN continuum. An 
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anti-correlation between the HI column density and the 
linear projected size of the radio source is now well 
established (jVermcul en et al. 2 003. Pihlstrom et al. 2003, 
|Gupta fc Saikia 2006[ ). Smaller sources (< 0.5 kpc) tend 
to have higher Hi column density than larger sources (> 
0.5 kpc). If not driven by uncertainties in the spin tem- 
perature of the HI absorbing gas, this anticorrelation can 
be explained by GPS galaxies hosting young radio sources, 
which evolve in a disk distribution of gas with a power-law 
radial density dependence. The same explanation could lie 
behind the (tentative) anti-correlation between X-ray col- 
umn density and radio size in GPS galaxies (Fig. [T0|). One 
can assume a scenario where the radio and X-ray source 
are seen through the same line-of-sight and are both em- 
bedded in a screen of ionised gas, responsible for X-ray 
photoelectric absorption as well as for radio free-free ab- 
sorption. From the definition of the ionisation parameter, 
it follows for the X-ray regime: 

where n e ^ is the plasma density in units of 10 3 cm -3 , Rk pc 
is the distance between the ionising source and the inner- 
most side of the absorbing material in kpc, and L44 is the 
ionising luminosity in units of 10 44 erg s _1 . It is plausi- 
ble that densities in GPS sources are high enough that 
free-free absorption may play a role (jO'Dea et al. 1991, 
IStawarz et al. 2008|) . The emission measure for free-free ab- 
sorption at 5 GHz is (see, e.g., Osterbrock 1977): 

n 2 e d pc ~ 90rT 4 135 

where d pc is the path length in parsecs and T4 is the temper- 
ature in units of 10 4 K. Assuming R ~ d, one can constrain 
the plasma temperature: 

T 4 <; (L 4 4?73o) 1 ' 5 -R A r p 2 c 2 

where 7730 is the conversion factor between the absorption- 
corrected 2-10 keV and the bolometric luminosity divided 
by 30. For standard values of ^30 (jElvis et al. 19 94b). al- 
most all GPS sources in our sample require temperatures 
above that at the dust sublimation radius in this scenario. 
This outcome could be tested by measurements of the op- 
tical reddening in our sample. 

If the compact jet "drills" its way through such a 
medium, would it be significantly decelerated or even "frus- 
trated" , i.e. permanently confined? This issue was discussed 
by Guainazzi et al. (2004) for the case of the Compton- 
thick absorber in OQ+208. They concluded that for large 
inclination angles or large thickness of the Compton-thick 
layer, the jet could have been significantly decelerated by 
the interaction with the ambient medium. Although perma- 
nent confinement is unlikely in OQ+208, underestimating 
the evolution time scale by one-two orders of magnitude 
is possible. We revise their argument for the whole sam- 
ple of GPS galaxies. The expansion time for a jet prop- 
agating in pressure equilibrium through in homogeneous 
medium can be expressed as (jScheuer 19741 [Carvalho 1985, 
IGarvalho 1998]) : 

t e ~ 6 rl 6 5 o pc L^% (cos z)" 1 - 5 ^ 5 ^ 5 s 

where Li„j^4 is the (unknown) luminosity injected in the 
jet in units of 10 44 erg s _1 , and fiio is the jet opening angle 



in units of 10 degrees. If we assume the censored best-fit 
for the linear size versus X-ray column density relation: 

t e ~ 4 x 10 5 r° 9 pc L-°j 5 44 years 

where we have compacted all the geometrical factors in the 
variable O = (cos i) -15 ^^ 5 - There are admittedly several 
uncertainties in deriving the above scaling law. Most im- 
portantly, we know very little of the actual gas distribution 
in the environment surrounding compact radio sources (as 
well as in any other type of AGN). Still, the above esti- 
mate indicates that a GPS may expand during a time of 
the order of at least 10 5 years. This is longer than currently 
existing observational estimates (Poladitis & Conway 2003; 
Gugliucci et al. 2005), but not enough to make the hypoth- 
esis of permanent confinement viable. 

6.5. Evolution of GPS sources 

If the jet can indeed survive its eventful youth, and grow to 
reach a full level of kpc-scale maturity and beyond, what 
would it look like? The radio-to-X-ray luminosity plane 
does not ultimately elucidate the possible connection be- 
tween GPS and "mature" radio galaxies. GPS galaxies are 
intriguingly well aligned along the extrapolation at high ra- 
dio power of the correlation between radio core and X-ray 
luminosity valid for FR I radio galaxies (Eva ns et al. 2 006) . 
This correlation was proposed as evidence supporting a jet 
origin of unabsorbed X-ray spectral components in FR I 
radio galaxies. On the other hand, GPS galaxies have a 
comparable X-ray luminosity to FR lis (cf. Fig. [TTj) . In 
FR lis, the obscured X-ray spectral component is probably 
associated with accretion onto the supermassive black hole 
obscured by a "torus" , in analogy to the accepted parad igm 
applicable to radio-quiet AGN (|Piconcelli et al. 2008[) . If 
the X-ray emission in GPS galaxies is due to accretion, the 
evolution of the radio and X-ray wavebands could be to- 
tally decoupled. Radio power would decline with the linear 
size (see, e.g., Fanti et al. 1995) while the sources expand 
through the ISM; at the same time the accretion disk would 
maintain a stable flow. At the end of their infancy, GPS 
galaxies would reach maturity as FR II radio galaxies. 

Evolutionary scenarios require that the radio luminosity 
of GPS sources decreases with evolution, not to exceed the 
number of observed FR II galaxies (Readhcad et al. 1996). 
Recently, Stawarz et al. (2008) have proposed an evolution- 
ary model for GPS sources, which explicitly predicts the 
dependency of the broadband Spectral Energy Distribution 
on the source linear size. In their model high-energy emis- 
sion is produced by upscattering of various photon fields by 
the lobes' electrons. They predict a decrease of the X-ray to 
radio luminosity ratio by 1-2 orders of magnitude when the 
GPS source size increases from 30 pc to 1 kpc. In Fig. [13] 
we compare this prediction with our observations. There is 
no evidence for a strong anti-correlation between the two 
quantities; the slope of the best linear fit is 0.3 ± 1.5 (Ter 
statistical error). Still, the observational data are consis- 
tent with the Stawarz et al. model predictions, if the GPS 
galaxies in our sample cover a wide range in jet kinetic 
power. 

It is still impossible with the available data to decide be- 
tween an evolutionary scenario in which the source evolves 
into a conventional FR II or FR I radio galaxy. To achieve 
this goal, photometric X-ray observations of sizable samples 
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Fig. 13. X-ray to radio luminosity ratio versus the linear 
size for the sources of our GPS sample. The lines indicate 
the predictions of the Stawarz et al. (2008) model for jet 
kinetic powers ranging from 10 44 erg s" 1 to 10 47 erg s _1 . 
Dashed lines: power-law injection function (Fig. 2 in their 
paper); dotted lines: broken power-law injection function 
(Fig. 3 in their paper). 
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would be crucial. Fortunately, it is now in principle possible 
to perform this experiment, thanks to the point-like source 
sensitivity and scheduling flexibility of Chandra. 

7. Summary and conclusions 

This paper reports our current knowledge on X-ray emis- 
sion in GPS galaxies. For the first time a complete radio- 
selected sample of GPS galaxies has been almost entirely 
observed with a modern X-ray observatory (mostly with 
XMM-Newton) . The sample is comprised all the z < 1 
sources of the Stanghellini et al. (1998) sample having a 
5 GHz flux density >1 Jy. The main results of our study 
can be summarised as follows: 

— we obtain a very large detection fraction; all the sources 
of our sample yield a detection in the soft X-ray band 
(0.5-2 keV), whereas 15 out of 16 are detected in the 
full band (0.5-10 keV) 

— in almost all simple power law modified by pho- 
toelectric absorption represents an adequate description 
of the 0.5-10 keV spectrum. In a few objects the number 
of net counts is not good enough to allow a full spectral 
analysis. In this case basic spectral parameters were de- 
rived from hardness ratios assuming the baseline model 
above 

— the mean of the spectral indices distribution is (r) = 
1.66 (err = 0.36). Although the uncertainty on this pa- 
rameter is still too large to pinpoint the physical mecha- 
nism responsible for the observed X-ray emission, at face 
value the T distribution is closer to those AGN classes 
whose X-ray emission is believed to be dominated by 



accretion: Seyfert Galaxies and the obscured spectral 
component in FR II radio galaxies 

— we report the possible discovery of a Compton- 
thick AGN in PKS1607+26. Radio-lou d Compton-thic k 
AGN are still an elusive population (Comastri 2004). 
Together with OQ+208 - the other Compton-thick AGN 
in the sample - PKS1607+26 is the only GPS galaxy 
where an X-ray emission line has been detected, possi- 
bly associated with Fe K Q fluorescence 

— X-ray spectra of GPS galaxies are significantly ob- 
scured. The mean value of the column density dis- 
tribution (without PKS1607+26, due to pending un- 
certainties on the identification of this source) is 
(Nft PS ) — 3 x 10 22 cm -2 with a standard devia- 
tion crjVff— 0.5 dex. Such a value is much larger than 
column densities measured in a control sample of 
FR I radio galaxies, but still less than column densi- 
ties covering accretion-related X-ray spectral compo- 
nents in FR II radio galaxies (fBalmavc rde et al. 20061 
IBelsole et al. 20061 lEvans et al."2006|) . 

— the X-ray column density measured in almost all GPS 
galaxies is larger than the HI column density mea- 
sured in the radio by a factor of 10 to 100. This 
could be the signature of phys i cally different absorp- 
tion systems ( Vi nk et al. 20051 IGua inazzi et al. 2006) 
or of a single system characterised by a spin temper- 
ature ^10 3 K (jOstorero et al. 2009|) . We report a pos- 
sible anti-correlation between the projected linear size 
of the radio source and the X-ray column density, anal- 
ogous to the anti-correlation between radio size and HI 
column density reported by Pihlstrom et al. (2003) and 
Gupta et al. (2006) 

— GPS galaxies occupy a specific locus in the radio versus 
X-ray luminosity plane. They lie well on the extrapo- 
lation to high radio powers of the correlation between 
these two quantities discovered in low-luminosity FR I 
radio galaxies. On the other hand, GPS galaxies exhibit 
a comparable X-ray luminosity to FR II radio galaxies, 
notwithstanding their much greater radio luminosity 

— GPS galaxies occupy the same locus as Seyfert galax- 
ies in the 0[III] to X-ray luminosity ratio versus X-ray 
column density diagnostic plane 

If the bulk of the X-ray emission in GPS galaxies is due 
to the accretion disk, one may expect that the X-ray and 
the radio source evolution in GPS galaxies would be de- 
coupled. The accretion disk would maintain a stable flow 
of gas and (mostly sublimated) dust to the supermassive 
black hole, while the radio source fades away with its expan- 
sion into the ISM. This would imply an increase of the X- 
ray to radio luminosity ratio with size. On the other hand, 
models of the dynamical evolution of GPS sources, where 
the X-ray emission is primarily produced by Compton up- 
scattering of ambient photons (Stawarz et al. 2008), pre- 
dict a decrease of the X-ray to radio luminosity ratio with 
size. The Lx / L^qhz ratios observed in our sample are con- 
sistent with both possibilities. Measurements in 7-rays by 
Fermi (as originally proposed bv lStawarz et al. 2008|) could 
be crucial to discriminate between them. Alternatively, the 
large scatter observed in this quantity may be indicative of 
a coupling between the accretion disk and the jet activity 
driven by disk instabilities. 

The evolutionary scenarios described above postulate 
that GPS sources are young objects, as also indicated by the 
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Nh° 


r n 




(10 20 cm" 2 ) 


(lO" 4 cm" 2 s" 1 ) 




2.29t£o6 7.40 ± 0.04 


1.00/350 



Spectral Energy Distribution 



a in excess of the Galactic contribution along the line-of-sight to 
XMMU J002200. 8+000655. 



direct measurement of their dynamical age. Eventually GPS 
galaxies would reach their full maturity as classical FR II 
radio galaxies. However, column densities <;10 22 cm~ 2 
fully surrounding the expanding radio source could sig- 
nificantly brake, if not entirely inhibit, this state, lead- 
ing to a significant underestimate of dynamical ages based 
on hotspot recession velocity measurements. Extending 
the number of X-ray measurements of low-luminosity 
{Lb GHz = 10 40 ~ 42 erg s _1 ) GPS galaxies is the next step 
we intend to pursue, in order to pinpoint the endpoint of 
their evolution. 



Appendix A: a serendipitous blazar in the field of 
4C+00.02 

A bright off-axis source is visible in the EPIC field 
of view of the XMM-Newton 4C+00.02 observation. 
This source will be referred to in the following as 
XMMU J002200.8+000655. It is outside the field of view 
of the XMM-Newton Optical Monitor. The best-fit param- 
eters when the baseline model is applied to its spectrum 
are summarised in Table The EW of a unresolved Fe 
K Q neutral fluorescent line is constrained to be lower than 
400 eV. 

This source was already known as 
1RXS J002200.9+000659 flVoges et al. 19991 ). We have 
searched with Aladin for available measurements at 
other wavelengths. We found da ta in GALEX GR4, SDSS 
(l Adelman-McCarthy et al. 2008[ ), 2MASS, and FIRST 
(White_elaIi99Z|) results. In Fig. [14] we compare the 
Spectral Energy Distribution (SED), with the average 
SED for radio-quiet and radio-loud quasars normalised 
at 1 keV (jElvis et al. 1994b|) and with the "blazar track" 
corresponding to objects with the same radio luminosity 
as XMMU J002200.8+000655 (|Fossati et al. 1998|) . As 
can be seen in FigQT] the blazar SED track qualitatively 
agrees with the XMMU J002200.8+000655 SED. This 
source appears in NED as a BL Lac candidate, and in 
Veron-Cetty & Veron (2006) as a confirmed BL Lac object. 
Its optical spectrum was previously studied by Collinge et 
al. (2005), but heavy contamination by the host galaxy 
prevented its properties from being properly studied. 



Appendix B: summary of the radio and X-ray prop- 
erties of the whole GPS sample 

In Table [6] we report a summary of the X-ray proper- 
ties of the whole GPS sample discussed in this paper, to- 
gether with the 5 GHz luminosity after O'Dea (1998) and 
Stanghellini et al. (1998). 
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Fig. 14. The Spectral Energy Distribution (SED) for 
XMMU J002200.8+000655. The photometric points are 
compared with standard SED for Blazars (boxes; 
IFossati et al. 1998)) and radio-quiet an d radio-loud quasar s 
(dotted line and solid line, respectively: lElvis et al. 1994a[) . 
normalised at the value of XMMU J002200.8+000655 at 
1 keV. 
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